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We present large Nf QCD/QED as a test bed for improved pressure calculations, 
show how to apply the hints obtained on optimized renormalization scales at large 
Nf to finite Nf = 2, and compare the results to recent lattice data. 



1. Introduction 

In the deconfined phase of QCD, strict perturbative calculations of ther- 
modynamic potentials show poor convergence when approaching the phase 
transition. There have been a number of attempts to overcome this prob- 
lem by reorganization or partial resummation of the perturbative expan- 
sions, like HTL perturbation theory^ or ^-derivable approximations for 2PI 
skeletons^, but so far independent verification of these models was only pos- 
sible through lattice simulations^. It is therefore instructive to consider the 
exactly solvable special case of large number of flavors (large Nf), in which 
the improvements above can be tested"*. The large Nf limit can further- 
more be easily extended to finite chemical potential^. Also, effects relevant 
to full QCD, like the anomalous specific heat at low temperatures, can be 
readily studied in the large- A^/ limit^'^. 



2. Large Nf 

The limit of large number of flavors is formed by sending Nf to infinity, 
while keeping the combination g^Nf as well as Nc of the order of 0(1). 
The diagrams contributing to the free energy are at leading order (LO) 
only the fermion loop, while at next-to-leading order (NLO) an infinite 
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Figure 1. Exact result for the interaction pressure at zero chemical potential as a func- 
tion of g'^ffiP-MS = ttT). The dashed line is the perturbative result, evaluated with 
renormalization scale pMS = Mfac ! the light band includes the numerically determined 
coefficient to order g^g (with its estimated error) also at ^fac- The result marked 
"sffi = deft" corresponds to choosing p,MS such that the order-g^g coefficient vanishes 
and retaining all higher-order terms contained in the plasmon term oc m^. 



number of ring diagrams consisting of a boson loop with any number of 
fermion loop insertions have to be summed up^. We can treat massless 
QCD and ultrarelativistic QED at the same time by defining an effective 
coupling as g^^ = for QCD and g^^ = e^Nf for QED. Large Nf 

contains a Landau pole of the order of Al ~ ^ exp(67r^/g^ff ), but the 
resulting ambiguity for the thermal pressure at NLO is suppressed by a 
factor (max(T, /i)/AL)'*. 

After subtracting off the vacuum part of the ring diagrams and applying 
Schwinger-Dyson resummation^, the NLO thermal pressure is given by 



^NLO f d^q 



Ng J (2^) 



[2 ([n, + i]Imln(g2 _ + Ht + H^ac) 

Jo 

-ilmln(g2 -q2 _^n^^c)) 



+ (k + i]Imln(i!^4t2|tIL:-) _ llmln( ^^;l+g- ))] (1) 

with the bosonic distribution function nb(aj) = l/{e^/'^ — 1) and the 
gauge-boson self energy functions IIt and 11^. These cannot be given 
in closed form except for their imaginary parts®, but are represented by 
one-dimensional integrals. We therefore have to evaluate the integrals nu- 
merically. Parts proportional to rif, can be safely integrated in Minkowski 



February 2, 2008 0:46 Proceedings Trim Size: 9in x 6in 



ipp 



3 

space, but terms without nt are potentially logarithmically divergent. We 
compute them by introducing a Euclidean invariant cutoff'*. 

3. Numerical results 

Figure 1 shows the numerical result for /i = as a function of .g^fj (/Ims = 
ttT). For small coupling the coefficients to order g^f^ which are not yet 
known analytically can be extracted numerically^. The large renormal- 
ization scale dependences of successive perturbative approximations to or- 
der g^g beyond g'^^ ^ 4 can be fixed by applying "fastest apparent con- 
vergence" (FAC) in the rn^ parameter of dimensional reduction. Using 
Mms = Mfac = 7re*/^~'''T we obtain good agreement up to g^^^ ^ 9. The 
result can be further improved by the procedures explained below Fig. 1. 




Figure 2. Exact result for the large- A^j interaction pressure P — Pq normalized to 
Ngiw'^T'^ + 1^1^)^ as a function of (/iivis) with /i^g = tt^T^ + and = arctan 

For non- vanishing chemical potential /i we use the fermionic distribution 
function 

which enters via the gauge boson self-energy expressions 11^ and 11^. In 
Fig. 2 we display our exact results for the interaction pressure P — Pq k 
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for the entire fi-T plane (but reasonably below the scale Landau pole) . The 
figure shows a kink at = 45° corresponding to fj, = ttT indicating that a 
simple scaling behavior for the pressure at small chemical potentials would 
break down at larger chemical potentials^. 




5 10 15 20 b i:j 15 20 



dad-') Kill/') 

Figure 3. Percentage errors of the perturbative result for the interaction part of the 
pressure to order g^f^ in the large-Af^ limit as a function ot = arctan(7rT//i) and 
9g{[{P-Ms) at pMS = TT-^T^ + /^^ for two choices of p,MS'- Fastest apparent convergence 
of P as well as m|, (FAC-m, left panel), and of g|, (FAC-g, right panel). The brightest 
area corresponds to an error of less than 1%, the darkest ones to an error of over 100%. 

Figure 3 shows a comparison between the perturbative result obtained 
by dimensional reduction^ through order with complete analytic depen- 
dence on arbitrary T and fi to the exact result in the large Nf limit. In 
the left panel the renormalization scale is /Ims = Mms~^~™- alterna- 
tive choice for the renormalization scale is to set /Ims such that the g'^g 
correction in the dimensional reduction coupling parameter g'^ is put to 
zero (FAC-g, right panel). The accuracy of the results is comparable in 
both cases and decreases slowly with increasing chemical potential, apart 
from an accidental zero of the error; for tp < 18°, i.e. T < 0.1/i, the errors 
eventually start to grow rapidly, marking the breakdown of dimensional 
reduction. This is precisely the region where non-Fermi-liquid effects lead 
to anomalous TlnT"-'^ terms in the entropy and specific heat^'^. 

In Fig. 4 we display the FAC optimized results^" at finite Nf — 2 for 
the difference AP — P{T,fj.) — P{T,0) for various /i/T corresponding to 
recent lattice resuhs" assuming Tq = T^=° = 0.49 Aqcd At T/Tq = 2 
our FAC-g and FAC-m results exceed the not-yet-continuum-extrapolated 
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Figure 4. The difference AP = P{T,fi) - P(r,0) divided by T"* for Nj =2 using the 
unexpanded three-loop result from dimensional reduction for fi/T = 0.2, . . . , 1.0 (bottom 
to top). Shaded areas correspond to a variation of fiyis around the FAC-m choice by a 
factor of 2; dashed and dash-dotted lines correspond to the FAC-g and FAC-m results, 
respectively. Also included are recent lattice data (not yet continuum-extrapolated!) 
assuming T^f^^ = 0.49 Aqcd- 



lattice data consistently by 10.5% and 9%, respectively, which is roughly 
the expected discretization error 
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